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ABSTRACT 

We suggest that polar rings and/or counter- rotating disks in flattened galaxies can be 
formed from stars captured at the Binney resonance, where the rate of precession of 
the angular momentum vector of a disk star equals the pattern speed of a triaxial halo. 
If the halo pattern speed is initially retrograde and slowly decays to zero, stars can be 
trapped as the Binney resonance sweeps past them, and levitated into polar orbits. If 
the halo pattern speed is initially retrograde and slowly changes to prograde, trapped 
stars can evolve from prograde to retrograde disk orbits. The stellar components of 
polar rings formed by this process should consist of two equal, counter-rotating star 
streams. 
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1 INTRODUCTION 



This paper is motivated by the remarkable discovery (Rubin, Graham & Kenney 1992) that the otherwise normal E7/S0 



galaxy NGC 4550 contains two stellar dis ks rotating in op posite directions. The two disks are similar in total luminosity and 
scale length and approximately coplanar ( Rix et al. 1992|) ; one is accompanied by an extended gas disk. 

A second case of a count er-rotating disk is the Sab g alaxy NGC 7217, in which 20-30% of the disk stars are on retrograde 
orbits, independent of radius ( Merrifield fc Kuijken 1994 ). Counter-rotating disks are rare: Kuijken, Fisher & Merrifield (1996) 
examined 28 SO galaxies and found no counter-rotating components containing more than ~ 5% of the total disk light. In 
contrast, roughly 20% of the gas disks in SO galaxies counter-rotate (Bertola, Buson & Zeilinger 1992); however, these gas 



disks are generally much smaller than the accompanying stellar disks, so even if they form stars (as in NGC 3593, Bertola et 
al. 1996) they are unlikely to produce two stellar disks of similar size as in NGC 4550 or NGC 7217. 

Several formation mechanisms for counter-rotating disks have been discussed, (i) A merger could add material on retro- 
grade orbits to a pre-existing stellar disk, but mergers are likely to overheat the original disk even in favourable cases where 
the merging galaxy is gas-rich or its orbit lies in the original disk plane (Phakar & Ryden 1998). (ii) Hierarchical models of 



galaxy formation predict that the mean angular momentum of infalling gas varies substantially during the lifetime of a galaxy, 
so that early infall could produce a gas disk that later forms stars, while late infall subsequently brings in retrograde gas 



that forms stars in turn. This proposal does not explain why the scale lengths of counter-rotating disks are similar (Thakar 
fe Ryden 19*9^ ), or why NGC 4550 satisfies the normal TuUy-Fisher relation ( [Rix et al. 1992| ). (iii) Evans & CoUett (1994) 
point out that box orbits in a triaxial potential can evolve into loop orbits if the potential slowly becomes more axisymmetric 
(an effect expected from late infall; Dubinski 1994). After this process stars will occupy both direct and retrograde loops — in 
precisely equal numbers if the triaxial potential is non-rotating — thereby naturally creating counter-rotating disksQ 

Polar-ring galaxies are early-type (usually SO) galaxies containing an outer ring of gas, dust and stars on orbits that are 
approximately circular and nearly perpendicular to the symmetry plane of the galaxy. At least 0.5% of SO galaxies have polar 



rings, although this is probably an underestimate because of orientation-dependent selection effects (Whitmore et al. 1990) 



The catalog by Whitmore et al. lists only six kinematically confirmed polar-ring galaxies but many more candidates without 
kinematic data. Like counter-rotating disks, polar rings are usually assumed to form from the merger of a companion galaxy 
or late gas infall (e.g. Katz & Rix 1992, Bekki 1998); if the gravitational potential of the primary galaxy is triaxial then there 



* The dynamical stability of counter-rotating disks is discussed by Sellwood &; Merritt (1994). 
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is a range of initial conditions from which dissipative material will settle into a polar orbit perpendicular to the long axis 



( steiman- Cameron & Durisen 1984; Thomas, Vine & Pearce 1994) 



In this paper we describe a novel way to form counter-rotating disks and/or polar rings. The inclinations of disk-star 
orbits can be excited by resonant coupling to a triaxial halo potential. The location of the relevant resonances is determined 
by the vertical and azimuthal frequencies, Q2 and (see and the pattern speed of the halo, Qp. In particular, Binney 
(1978, 1981) has stressed the importance of the resonance at 

Q3-n2= fip, (1) 

which we call the Binney resonance. The Binney resonance occurs when the pattern speed matches the rate of precession of 
the angular momentum vector or node, tj = SI3 — f22. For low-inclination orbits in oblate potentials, ^2 > ^3, so the Binney 
resonance occurs for retrograde pattern speed. 

We suppose that the disk is embedded in a triaxial halo that initially rotates with a retrograde pattern speed i}pi < 0. 
The pattern speed is expected to change slowly as the halo acquires dark matter by infall; we assume that the pattern speed 
increases, reaching a final value Qpf > 0. If f2p changes sufficiently slowly, stars with small initial inclinations io are trapped in 



the Binney resonance as f2p sweeps past 573 — 02. As f2p increases further, the orbits of the trapped stars are levitated (Sridhar 



& Tot ma 1996) to higher inclination while remaining nearly circular, becoming polar orbits as Qp crosses zero. Thus if \Qp\ 
gradually decays to zero, an outer polar ring is formed from disk stars with 572 — fis < \i^pi\- If ^p crosses zero to positive 
values, the trapped stars evolve onto retrograde orbits, until they are eventually released from the Binney resonance when Qp 
sweeps through —Q2 + fl3>0. The inclinations after release are near tt — io, so the fiipped stars form a counter- rotating disk 
with the same thickness as the initial disk. 

A close cousin of this process was already discussed by Heisler, Merritt & Schwarzschild (1982; see also van Albada, 
Kotanyi & Schwarzschild 1982). They recognized that a sequence of "anomalous" inclined orbits bifurcated from the closed 
short-axis loop orbits at the Binney resonance in a rotating triaxial potential, and speculated that gas might evolve onto 
the anomalous orbits as a result of dissipation. However, they focused on a sequence of orbits of decreasing energy at fixed 
pattern speed, which terminates in a short-axis orbit that cannot be occupied by collisional material. In contrast, we examine 
the behavior of dissipationless material in a system with time-varying pattern speed. Our mechanism is also related to the 
levitation process discussed by Sridhar & Touma (1996), who focus on a different resonance as a means of forming a thick 
disk. 

In §^ we describe the classification of resonances in nearly axisymmetric potentials and justify our focus on the Binney 
resonance. In §^ we investigate capture and release of stars at the Binney resonance using a simplified dynamical model. §^ 
describes the results of numerical integrations, and ^ contains a discussion. 



2 RESONANCE CLASSIFICATION 

To determine the most promising sites for resonant capture, we first consider integrable motion in an axisymmetric potential. 
We define action-angle variables (I, w), such that Hq = -ffo(I) is the axisymmetric Hamiltonian and the trajectory of a particle 
is given by 

I — const, w = const -I- fit, (2) 
where 

^^(I) - (3) 

The actions in a spherical potential can be chosen as follows: Ji is the radial action, which is zero for circular orbits; I2 is 
the vertical or latitudinal action, which is zero for prograde equatorial orbits and in general is equal to J — Jz, where J is the 
total angular momentum and Jz its z-component; I3 is the azimuthal action, which is equal to Jz- Note that 7i > 0, I2 > 0, 
/3 < hh. 



The geometrical interpretation of these actions remains similar for integrable axisymmetric Hamiltonians ( de Zeeuw 



1985 ); in particular, la is still equal to Jz, prograde equatorial orbits still have I2 = 0, and the analogs of circular orbits 
(Ji = 0) are shell orbits, which occupy a two-dimensional axisymmetric surface of zero thickness. 
Now we add a weak, non-axisymmetric perturbing potential of the form 

Uiiv,t)=A{T,e)cos{m[<j,^Mm, (4) 

where m 7^ is a positive integer and (pp = Sip is the pattern speed. In action-angle variables this potential can be written as 
(e.g. Tremaine & Weinberg 1984) 

Ui(l,w,t) = ^Ai^(l)cos[k-w~m(Pp{t)] (5) 

k 
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where k is an integer triple, with fcs = m. Terms in the perturbing potential (j^ with — conserve the azimuthal action 
^3 = Jz', terms with fei = conserve the radial action /i; and terms with ^2 = conserve the vertical action I2 ~ J — Jz- 
When 7i = the radial angle wi is undefined, so ylk(/i = 0) must vanish for fci 7^ 0; more generally ylk ^ /l*^!!/^ as Ji 0. 
Similarly — > Ij*^'^'^ as 72 ^ (in celestial mechanics these constraints are called the d'Alembert characteristic). We shall 
assume that the perturbing potential is symmetric around the equatorial plane, which requires that j4k ~ unless ^2 is even. 

We are interested in the resonant excitation of inclination in low-inclination, near-circular disk orbits. Resonance occurs 
when the rate of change of the resonant angle k ■ w — m(j)p{t) is zero, yielding the resonance condition 

k ■ f2 = mfip. (6) 

Since Ak ^ as Ji ^ when ki 7^ 0, the strongest inclination excitation for near-circular orbits comes from terms with ki = 0. 
We shall also focus on terms with k^ = m — 2, since barlike and triaxial perturbations are the strongest non-axisymmetric 
features in galaxies. Since Ak tx near-equatorial orbits are most strongly affected by terms with k2 = 0, ±2 (recall 

that k2 must be even); however, resonances with k2 — cannot excite inclination (for these resonances W2 is ignorable, so the 
vertical action 72 is conserved), so we restrict our attention to k2 = ±2. For the triples k = (0, ±2, 2) the resonance condition 
(|) is 

±n2 + i^3 = ^p- (7) 

In potentials that are not highly flattened, the vertical and azimuthal frequencies are similar, ^2 ~ ^^3- If the pattern speed is 
slow, \Q,p\ < Sl2,3 (as would be expected for the perturbation from a triaxial halo) resonance is therefore more likely to occur 
for k2 = —2 than for k2 = 2. Thus we are led naturally to the Binney resonance k = (0, —2, 2) (eq. |l|). 

At the Binney resonance 7i is adiabatically invariant because the perturbations with fci 7^ are rapidly varying. Moreover 
W2 and 103 appear in the potential only in the combination 2{ws — W2), so 72 -1- 73 is conserved. In nearly spherical potentials 
72 -1- 73 ~ J, so the total angular momentum is approximately conserved. In other words the resonance affects only the 
inclination of the orbit, and not its total angular momentum or eccentricity. In effect, the evolution of a circular orbit near 
the Binney resonance is that of a rigid spinning hoop with the same radius and angular momentum, an analogy that we shall 
pursue in the next section. 



3 A SIMPLIFIED MODEL FOR RESONANCE CAPTURE 

We consider motion in a rotating, triaxial potential of the form 

f/(r,t) = $(^fVg + ^J); (8) 
here the rotating coordinates {x,y) are related to the inertial coordinates {x,y) by 

x + iy={x + iy)exp{i(j}p), '^=Q.p{t), (9) 

where Q.p is the pattern speed. We assume that the a;-axis is the long axis, so < p, g < 1, and that the symmetry plane of 
the disk is the x-y plane. Since loop orbits around the intermediate axis are unstable, the 2-axis must be the short axis, so 
q <p. Thus 

0<q<p<l. (10) 

We assume that the axis ratios are not too far from unity — typically we choose q = 0.8, p = 0.9 — so we can expand the 
potential as 

C/(r, t) = $(r2) + ^ [(p-2 - l)y^ + {q-^ - 1)^^] ; (11) 

here we have replaced d^/dr^ by ^V^/r^ where V{r) is the circular speed at radius r. 
For simplicity we assume that the time variation of the pattern speed is given by 

ftp{t) = Clpt, (12) 

where Clp is a constant. 

Since we are interested in nearly circular orbits, we replace the particle by a circular hoop of radius r; the angular mo- 
mentum of the hoop is a constant J and its orientation is specified by the canonical momentum-coordinate pair (Jz, tj), where 
Jz = J cos i is the z-component of the angular momentum, i is the inclination, and uj is the longitude of the ascending node. 
The Hamiltonian of the hoop is H{Jz,u!,t) = {U{r,t)), where {•) denotes an average over the hoop. Neglecting unimportant 
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constants we have 
1 

"4 



H{Jz,ijJ,t) = — — l)[sin^((jj — (;/>p) + cos^ i cos^(cj — (?!>p)] + — 1) sin^ i} 



' 8 



(13) 



Throughout this analysis we can consider V to be constant, since the radius of the hoop is fixed. 

We now convert to new canonical variables {W,w) = (Jz/J, a; — (j>p) and a dimensionless time s — t/{Jl3), where /3 is a 
constant to be chosen below. Using equation (^2|) for the time-dependence of the pattern speed, the new Hamiltonian is found 
to be 

H{W,w,s) = _^(2g-2-p-2_i)M/2 

-^(p^^ - 1)(1 - W^)cos2w - {f3^J^np)Ws. (14) 
8 

We choose the timescale parameter /3 so that the coefficient of is — ^ ; thus 



2 ' 

H{W,w,s) = -^W^ -a{l - W^)cos2w~-yWs, (15) 
where 



/3 = 



2(2g-2 - 1)' 

4 

V2(2g-2 -p-2 - 1) 
16r2f2„ 

^2' 



(16) 



' l/2(2q-2 -p-2 - 1)2 

here we have replaced J by rl/, its value when p and q are near unity. The constraints (^o[) imply that 

0<a<i, /3>0, (17) 

and sgn(7) = sgn(f2p); for the axis ratios q = 0.8, p = 0.9, we have a = 0.1317. 

The Hamiltonian (^|) depends on two parameters: a, which determines the strength of the non-axisymmetric potential, 
and 7, which determines the speed of the resonance passage. In the limit where the potential is nearly spherical, 1 — p, 1 — q — 
0(e) ^ 1, we have a = 0(1), (3 = 0{e-^), 7 = 0(f2pe-^). The phase space is a sphere of unit radius with a simple physical 
interpretation: the angular-momentum vector J has azimuth w + <j>p — ^vr and colatitude i = cos~^ W . 

We are interested in the case where the pattern speed changes slowly, \Clp\ <C or |7j <C 1. In this case, over short times 
the trajectory closely follows a level curve of the Hamiltonian for fixed s, H(W,w, s) = h{s) (the "guiding trajectory"). The 
nature of the guiding trajectories depends on the topology of the level surfaces of H. In describing this topology we restrict 
ourselves to < a < ^, as required by equation (^); we may then distinguish the following stages (Figure ^: 

(a) 7s < — 2q — 1: There are stable equilibrium points at W — 1 (north pole) and W — —1 (south pole); the angle w 
circulates for all orbits. 

(b) —2a! — 1 < 7s < 2a — 1: There is a stable equilibrium point at the south pole, and an unstable equilibrium or saddle 
point at the north pole. There are also stable equihbria at w — i^vr, W = — 7s/(l-|-2Qf). The separatrix orbit passing through 
the north pole has energy hscp = — \ — is, and divides circulating orbits with h < hsap from librating orbits with h > hscp- 

(c) 2a — 1 < 7s < —2a + 1: There are stable equilibrium points at both poles, as well as stable equilibria at w — i^vr, 
W = — 7s/(l + 2a). In addition there are unstable equilibria at 10 = {0, vr}, W = — 7s/(l — 2a). The separatrix orbit passing 
through the unstable equilibria has energy hscp ~ |(7s)^/(l — 2a) — a, and divides circulating orbits with h < hscp from 
librating orbits with h > hscp- 

(d) —2a -|- 1 < 7s < 2a -I- 1: This is identical to stage (b) after the transformation 7s — > —7s, W — > —W. There are stable 
and unstable equilibria at the north and south poles respectively, and stable equilibria at w = if"", W — —75/(1 -I- 2a). 

(e) 2a -|- 1 < 7s: This is identical to stage (a) after the transformation 7s — > —7s, W — > —W: there are stable equilibria at 
W — H and w circulates for all orbits. 

Adiabatic invariance ensures that over long times the guiding trajectory evolves so as to preserve the action, which is 
(27r)~^ times the area on the phase-space sphere enclosed by the guiding trajectory (e.g Peale 1976, Henrard 1982, Borderies 
& Goldreich 1984, Engels & Henrard 1994). 

Let us follow the evolution of the orbit for the case 7 > 0, corresponding to a pattern speed that is initially negative 
(retrograde) but increasing. At large negative time, 7s <C —1, the Hamiltonian ( |l5| ) is dominated by the term —yWs and the 
guiding trajectory is W = Wo =constant, where Wo ~ cosio and io is the initial inclination. The initial action is (27r)~^ times 
the area of the north polar cap on the phase-space sphere that is enclosed by the initial trajectory, and is equal to 1 — Wo- 
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(a) a = 0.3, 7s = -4.0 




(b) a = 0.3, 7S = -0.7 




(c) a = 0.3, 7s = -0.2 



Figure 1. Topology of the level surfaces of the Hamiltonian (|l5|), for o < i. The plots show equal-area Aitoff projections of the sphere 
with longitude w and colatitude i = cos~^ W, where i is the inclination and w G [— 7r,7r] is the longitude of the ascending node in the 
frame rotating at the pattern speed Op. The plots show only stages (a)-(c); stages (d) and (e) are upside-down versions of (b) and (a) 
respectively. 



As the time s increases, the topology of the Hamiltonian eventually changes from stage (a) to stage (b). In the initial 
phases of stage (b), the guiding trajectory continues to circulate, despite the growing libration zones defined by the separatrix 
orbit through the north pole. Eventually the area occupied by the libration zones grows to 2n{l — Wo), so the action can 
no longer be conserved if the orbit continues to circulate. At this point the guiding trajectory crosses the separatrix and is 
captured into one of the two libration zones. 

The area occupied by the libration zones in stage (b) can be evaluated analytically, 

Aijs)^- '^'^1.,,., COS'' l+7g-4«'' ^4eQg-i -1-7S -2q - 1< 7s < 2q - 1. (18) 

^' ^ (l-4a2)i/2 2a7s 2a ' - / - v ; 

As the time s increases, A{'ys) increases from zero at the onset of stage (b) to Amax at the onset of stage (c), where 

•l-2a\i/2l 



^max = A{2a - 1) = 47r 



1 - 



/l-2ay 
Vl + 2q/ 



(19) 



Capture into libration occurs when A = 27r(l — Wo) and is certain if Amax > 27r(l — Wo); in other words capture is certain if 
'l-2a\i/2 



cosio > 2 



1. 



a + 2a , 
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Figure 2. Initial inclination iq versus final inclination for initially circular orbits with radius r = 1 in the triaxial potential defined by 
equation (p^. The evolution of the pattern speed is given by equation (^^, with Qpf = ~flp 
inclination is ~ 180° — in. 



: 0.4 and r = SOOO.For io ^ 50°, the final 



For our nominal value a — 0.1317, capture is certain if the initial inclination io < 58.2°. For larger inclinations, capture is 
probabilistic because the guiding trajectory crosses the separatrix orbit during stage (c), where transition to either libration 
or circulation can occur. In this case the capture probability can be computed using methods described by Henrard (1982). 

The captured orbits remain librating through stage (c) and in the initial phases of stage (d). Eventually they re-cross 
the separatrix into circulation; by symmetry their final inclination is just if — tt ~ io. In other words, a disk of stars in direct 
rotation is fiipped into a disk with the same radial profile and thickness but in retrograde rotation. This mechanism operates 
if (i) The initial pattern speed i}pi corresponds to stage (a) and the final pattern speed Sip/ to stage (e); this requires 



V . _2 

2^(^ 



1), 



V 

fip/ > — (q 



!)• 



(21) 



This condition can be restated in terms of uj — Qs — Q2, the nodal precession rate for low- inclination orbits, as flpi < cu, 
Qpf > —uj. (ii) The pattern speed changes slowly enough that the action is adiabatically invariant except near the separatrix. 
(iii) The initial inclination of the disk stars is sufficiently small (eq. po[ ). 

Another interesting outcome occurs if the pattern speed is initially retrograde and slowly decays to zero. In this case the 
stars will be trapped in polar orbits (stage (c) with s = 0); the trapping process populates the two separatrices equally, so 
the resulting polar ring will itself form a counter-rotating disk, perpendicular to the long axis of the triaxial potential. 



4 NUMERICAL RESULTS 

We have followed the evolution of test-particle orbits in a rotating triaxial potential of the form 



U{v,t) = IV' log + ^ + 



-2 , f . z' 



(22) 



where y is a constant and x and y are defined by equation (|9|). We take axis ratios q = 0.8, p = 0.9 and set V = 1. The 
particle is initially on a circular orbit with radius r = 1. The pattern speed is assumed to vary as 



Qpit) = fip, + {^p, - Qpf) (e-'/" - 1) 



(23) 



and the orbits are followed from t — to t — 4r. With these parameters the star can undergo resonant capture and release if 
(eq. §) 



Qpi < -0.281, Qpf > 0.281. 



(24) 
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Figure 3. Final energy and angular momentum for the flipped particles in Figure y (iq < 40°). The initial energy and angular momentum 
are shown by dashed lines. The changes in both quantities are small, confirming that the size and shape of the orbit are not changed by 
the Binney resonance. 



To illustrate resonant capture and release in the adiabatic limit, we have numerically integrated the orbits of 200 test 
particles with random orbital phases and nodes, varying the pattern speed according to the parameters ilpi = —0.4, Qpf — 
+0.4, T = 5000. The results are shown in Figure]^: stars with initial inclination io ^ 50° are flipped to inclination ~ 180° —io, 
while for io ^ 50-60° the final inclination exhibits a wide spread. These results are consistent with the conclusion from §^that 
resonant capture was certain for these parameters if io < 58.2°. Figure ^ shows that the fractional changes in total angular 
momentum and energy of the flipped particles are 0.1, confirming that capture in the Binney resonance changes only the 
inclination of the particle orbits, not their size or shape. 

To investigate the validity of the adiabatic approximation, we have integrated 200 test particles with r varying from 
100 to 2000, and initial inclinations distributed as n{io)dio oc iodio exp( — ^io/cr^), with a = 5°. The distribution of final 
inclinations shown in Figure ^ shows that most stars are captured when r <; 300; in physical units this corresponds to 
r ^ 7.5Gyr(r/5kpc)(200kms-VV). 

We have also investigated the case where the pattern speed is initially retrograde and slowly decays to zero. Figure ^ 
shows the distribution of final inclinations as a function of the timescale t, for Qpi = —0.4, flpf = 0. The initial inclination 
range was io — 0-10°. For r ^ 100 most of the particles are trapped in orbits with inclination near 90°, thus forming a polar 
ring. The angular-momentum vectors are aligned with the long axis of the potential, with equal numbers of stars having 
Jx > and Jx < 0. Thus the polar ring is itself a counter-rotating disk. 

Finally, we have examined whether resonance capture can occur in disk galaxies, by integrating test-particle orbits in a 
non-axisymmetric Miyamoto-Nagai potential 

GM 



<i>{R,z) = 



(25) 



(i2 + y2/p2 4. (a + (f,2 + ^2)1/2)2)1/2 ' 

where x and y are defined in equation (^) and a — 0.5, b — 0.1. The particles are initially in circular orbits of unit radius, 
and the pattern speed varies according to (|2^) with Qpi — —2, Qpf = 2. We find that for initial inclination io ^ 6°, most 
particles are captured and flipped into retrograde orbits with inclination 180° — io- The maximum inclination of the flipped 
orbits depends on the parameter 6, which controls the thickness of the disk mass distribution; for b = 0.2 most particles are 
flipped if io ^ 15-20°. 



5 DISCUSSION 

In rotating triaxial potentials with time- varying pattern speed, stars on short-axis loop orbits can be captured at the Binney 
resonance, where the pattern speed matches the precession rate uj of the angular momentum vector. As the pattern speed 
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Fig ure 4. (top) Final inclination for initially circular orbits with radius r = 1 in the triaxial potential defined by equations (E3) and 
(bsp, with = —0.4, f2pj- = 0.4, and t varying from 100 to 2000. The initial inclinations are distributed as dN oc jq exp(— iiy/a^)djoi 
where cr = 5° . The final inclinations are clustered near 180° for T > 300. 
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Figure 5. (top) Final inclination for initially circular orbits with radius r = 1 in the triaxial potential defined by equations (g^ and 
(bsh, with Q,pi = —0.4, Op J- = 0, and t varying from 10 to 10^. The final inclinations are clustered near 90°, so the captured stars form 
a polar ring, (bottom) Distribution of the azimuthal angle of the final angular momentum vector in the rotating frame. The angular 
momentum vectors are centered on the long or x-axis of the potential, with equal numbers of stars rotating in opposite directions. 
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Figure 6. Initial inclination ig versus final inclination for initially circular orbits with radius r = 1 in the Miyamoto-Nagai potential ( |25| ) 
with a = 0.5, b = 0.1. The evolution of the pattern speed is given by equation (p^), with Qpf = — flpi = 2 and t = 2000. For io 5i 6°, 
the final orbits are retrograde and equatorial. 



continues to change, captured stars will be carried to high inclinations without large changes in energy or total angular 
momentum. 

Binney's (1981) original discussion focused on the linear response of stars close to the resonance when the pattern speed 
was fixed. In contrast, by exploiting the theory of slow (adiabatic) resonant capture (Peale 1976; Henrard 1982; Borderies & 
Goldrych 1984; Engels & Henrard 1994), we can follow the nonlinear orbital evolution so long as the variation in pattern 
speed is slow enough. 

It is striking that slow resonant capture and subsequent inclination growth is certain for stars with inclination i ~ 0, no 
matter how weak the triaxial potential may be, even though i = is a formal solution of the equations of motion for the model 
Hamiltonian (^3|) for all time. The resolution of this apparent paradox is that (i) a Mathieu-type linear inclination instability 
is always present when the star is sufficiently close to resonance (Binney 1981); (ii) as the strength of the non-axisymmetry 
approaches zero, the drift rate of the pattern speed must also approach zero in order that the adiabatic approximation is 
valid. 

If the final pattern speed of the triaxial potential is near zero, stars captured into the Binney resonance will form a polar 
ring of long-axis loop orbits. This model explains naturally why polar rings do not extend to the center of the galaxy: only 
stars whose initial precession rate |Lj| is less than the initial pattern speed \Qpi\ can be captured. The stellar component of 
polar rings formed by this mechanism should exhibit two equal, counter-rotating star stream, a testable prediction. 

If the final pattern speed is positive and larger than the precession rate of the angular momentum vector, Qpf > \uj\, 
stars captured into the resonance with initial inclination io will be flipped onto retrograde orbits and then released with final 
inclination n — io. In disk potentials only the low- inclination fraction of the disk orbits are captured and flipped (Fig. so 
both the direct and retrograde disk can be composed of pre-existing stars. For mildly triaxial potentials, capture is certain 
for orbits with small or moderate inclinations; in this case we must rely on subsequent star formation to re-form the parent 
disk. This model explains naturally why the two components of the counter-rotating disks in NGC 4550 and NGC 7217 have 
similar scale lengths. It does not explain why the two components in NGC 4550 have similar luminosity but then the two 
components in NGC 7217 do not, and in any case there are strong selection effects that favor the discovery of counter-rotating 
disks with similar luminosity. 

A concern with this model is whether the required timescale for variation of the pattern speed is unrealistically slow. 
The time unit in our simulations is 2.4 x 10^ yr(r/5 kpc)(200 kms~^/l/), so timescales r ^ 400 may exceed 10^" yr, the 
natural timescale for variations in halo pattern speed. This should be compared with r ^ 300 required to flip orbits in the 
triaxial potential (Fig. ^ and r <; 1000 to flip orbits in the Miyamoto-Nagai disk (Fig. On the other hand we have used 
over-simplified model potentials and have not explored parameter space systematically. 
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Another issue is whether gravitational noise due to molecular clouds or spiral arms degrades the effectiveness of resonant 
capture, although this is not a problem for the majority of polar rings that are found in SO galaxies. 

An unresolved question is whether significant quantities of gas can be captured into the Binney resonance. The orbits 
of particles in the two libration zones around w — ii^-K intersect, so gas clouds in the two zones will collide at high speed, 
leading to rapid energy dissipation. There are thus two possibilities: either no gas is captured, or all of the gas is captured 
into one of the libration zones, leaving the other vacant. Numerical simulations are the best way to determine which of these 
two outcomes is more realistic]^. 

This research was supported in part by NSF Grant AST-9900316 and NASA Grant NAG5-7066. 
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t Numerica l simulations of polar rings (Habe & Ikeuchi 1985; Varnas 1990; Quinn 199l[ |Katz fc Rix 199^ ; [christodoulou ct al. 1992 



Bekki 1998) show that gas rings on the long-axis loop orbits in triaxial potentials can be stable in the presence of dissipation, even 



though the libration zone surrounds a local maximum of the averaged Hamiltonian (|13|). 
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